Molecular Simulation of Argon Adsorption and Diffusion in a Microporous Carbon with Poroelastic Couplings.

Neglected for a long time in molecular simulations of fluid adsorption and transport in microporous carbons, adsorption-induced deformations of the matrix have recently been shown to have important effects on both sorption isotherms and diffusion coefficients. Here we investigate in detail the behavior of a recently proposed 3D-connected mature kerogen model, as a generic model of aromatic microporous carbon with atomic H/C ∼ 0.5, in both chemical and mechanical equilibrium with argon at 243 K over an extended pressure range. We show that under these conditions the material exhibits some viscoelasticity, and simulations of hundreds of nanoseconds are required to accurately determine the equilibrium volumes and sorption loadings. We also show that neglecting matrix internal deformations and swelling can lead to underestimations of the loading by up to 19% (swelling only) and 28% (swelling and internal deformations). The volume of the matrix is shown to increase up to about 8% at the largest pressure considered (210 MPa), which induces an increase of about 33% of both pore volume and specific surface area via the creation of additional pores, yet does not significantly change the normalized pore size distribution. Volume swelling is also rationalized by using a well-known linearized microporomechanical model. Finally, we show that self-diffusivity decreases with applied pressure, following an almost perfectly linear evolution with the free volume. Quantitatively, neglecting swelling and internal deformations tends to reduce the computed self-diffusivities.

Simple and efficient algorithms based on Volterra equations to compute memory kernels and projected cross-correlation functions from molecular dynamics.

Starting from the orthogonal dynamics of any given set of variables with respect to the projection variable used to derive the Mori-Zwanzig equation, a set of coupled Volterra equations is obtained that relate the projected time correlation functions between all the variables of interest. This set of equations can be solved using standard numerical inversion methods for Volterra equations, leading to a very convenient yet efficient strategy to obtain any projected time correlation function or contribution to the memory kernel entering a generalized Langevin equation. Using this strategy, the memory kernel related to the diffusion of tagged particles in a bulk Lennard-Jones fluid is investigated up to the long-term regime to show that the repulsive-attractive cross-contribution to memory effects represents a small but non-zero contribution to the self-diffusion coefficient.

Free Volume Model for Transport in Flexible Kerogen of Source Rock's Organic Matter.

We build a model of transport of adsorbed fluid within the microporous network of kerogen's porosity, especially accounting for the adsorption-induced swelling exhibited by flexible kerogen structures. This model, based on Fujita-Kishimoto free volume theory that was historically developed for swellable polymers, is built over extensive results for the self-diffusion coefficients obtained by molecular dynamics calculations for a representative molecular model of kerogen designed to study the importance of flexibility effects on the properties of kerogen. To do so, we first highlight that transport within flexible kerogen incorporating the coupling between the dynamics of the fluid molecules and the kerogen matrix atoms does not introduce any significant collective effects in the usual long time limit. Then, we show that despite the slightly anisotropic diffusion properties, averaging over all the dimensions can still be performed in order to model the behavior of the transport properties with the amount of adsorbed fluid. Lastly, we link the increase of the self-diffusion coefficients and that of the accessible free volume with the fluid loading via the Fujita-Kishimoto model. We conclude by commenting on the evolution and significance of the model parameters over a broad range of thermophysical conditions.

Nonlocal Screening Dictates the Radiative Lifetimes of Excitations in Lead Halide Perovskites.

We use path integral molecular dynamics simulations and theory to elucidate the interactions between charge carriers, as mediated by a lead halide perovskite lattice. We find that the charge-lattice coupling of MAPbI3 results in a repulsive interaction between electrons and holes at intermediate distances. The effective interaction is understood using a Gaussian field theory, whereby the underlying soft, polar lattice contributes a nonlocal screening between quasiparticles. Path integral calculations of this nonlocal screening model are used to rationalize the small exciton binding energy and low radiative recombination rate observed experimentally and are compared to traditional Wannier-Mott and Fröhlich models, which fail to do so. These results clarify the origin of the high power conversion efficiencies in lead halide perovskites. Emergent repulsive electron-hole interactions provide a design principle for optimizing soft, polar semiconductors.

Solid-State Ionic Rectification in Perovskite Nanowire Heterostructures.

Halide perovskites have attracted increasing research attention with regard to their potential for optoelectronic applications. Because of its low activation energy, ion migration is implicated in the long-term stability and many unusual transport behaviors of halide perovskite devices. However, direct observation and precise control of the ionic transport in halide perovskite crystals remain challenging. Here, we have designed an axial CsPbBr3-CsPbCl3 nanowire heterostructure, in which electric-field-induced halide ion migration was clearly visualized and quantified. We demonstrated that halide ion migration is dependent on the applied electric field and exhibits ionic rectification in this solid-state system, which is due to the nonuniform distribution of the ionic vacancies in the nanowire that results from a competition between electrical screening and their creation/destruction at the electrodes' interfaces. The asymmetric heterostructure characteristics add an additional knob to control the ion movement in the design of advanced ionic circuits with halide perovskites as building blocks.

Methane Diffusion in a Flexible Kerogen Matrix.

It has been recognized that the microporosity of shale organic matter, especially that of kerogen, strongly affects the hydrocarbon recovery process from unconventional reservoirs. So far, the numerical studies on hydrocarbon transport through the microporous phase of kerogen have neglected the effect of poromechanics, that is, the adsorption-induced deformations, by considering kerogen as a frozen, nondeformable, matrix. Here, we use molecular dynamics simulations to investigate methane diffusion in an immature (i.e., with high H/C ratio) kerogen matrix, while explicitly accounting for adsorption-induced swelling and internal matricial motions, covering both phonons and nonperiodic internal deformations. However, in the usual frozen matrix approximation, diffusivity decreases with increasing fluid loading, as evidenced by a loss of free volume, accounting for adsorption-induced swelling that gives rise to an increase in free volume and, hence, in diffusivity. The obtained trend is further rationalized using a Fujita-Kishimoto free volume theory initially developed in the context of diffusion in swelling polymers. We also quantify the enhancing effect of the matrix internal motions (i.e., at fixed volume) and show that it roughly gives an order of magnitude increase in diffusivity with respect to a frozen matrix, thanks to fluctuations in the pore connectivity. We eventually discuss the possible implications of this work to explain the productivity slowdown of hydrocarbon recovery from shale immature reservoirs.

Intrinsic anion diffusivity in lead halide perovskites is facilitated by a soft lattice.

Facile ionic transport in lead halide perovskites plays a critical role in device performance. Understanding the microscopic origins of high ionic conductivities has been complicated by indirect measurements and sample microstructural heterogeneities. Here, we report the direct visualization of halide anion interdiffusion in CsPbCl3-CsPbBr3 single crystalline perovskite nanowire heterojunctions using wide-field and confocal photoluminescence measurements. The combination of nanoscale imaging techniques with these single crystalline materials allows us to measure intrinsic anionic lattice diffusivities, free from complications of microscale inhomogeneity. Halide diffusivities were found to be between 10-13 and ∼10-12 cm2/second at about 100 °C, which are several orders of magnitudes lower than those reported in polycrystalline thin films. Spatially resolved photoluminescence lifetimes and surface potential measurements provide evidence of the central role of halide vacancies in facilitating ionic diffusion. Vacancy formation free energies computed from molecular simulation are small due to the easily deformable perovskite lattice, accounting for the high equilibrium vacancy concentration. Furthermore, molecular simulations suggest that ionic motion is facilitated by low-frequency lattice modes, resulting in low activation barriers for vacancy-mediated transport. This work elucidates the intrinsic solid-state ion diffusion mechanisms in this class of semisoft materials and offers guidelines for engineering materials with long-term stability in functional devices.

Mesoscale structure, mechanics, and transport properties of source rocks' organic pore networks.

Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores' topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.

Poroelasticity of Methane-Loaded Mature and Immature Kerogen from Molecular Simulations.

While hydrocarbon expulsion from kerogen is certainly the key step in shale oil/gas recovery, the poromechanical couplings governing this desorption process, taking place under a significant pressure gradient, are still poorly understood. Especially, most molecular simulation investigations of hydrocarbon adsorption and transport in kerogen have so far been performed under the rigid matrix approximation, implying that the pore space is independent of pressure, temperature, and fluid loading, or in other words, neglecting poromechanics. Here, using two hydrogenated porous carbon models as proxies for immature and overmature kerogen, that is, highly aliphatic hydrogen-rich vs highly aromatic hydrogen-poor models, we perform an extensive molecular-dynamics-based investigation of the evolution of the poroelastic properties of those matrices with respect to temperature, external pressure, and methane loading as a prototype alkane molecule. The rigid matrix approximation is shown to hold reasonably well for overmature kerogen even though accounting for flexibility has allowed us to observe the well-known small volume contraction at low fluid loading and temperature. Our results demonstrate that immature kerogen is highly deformable. Within the ranges of conditions considered in this work, its density can double and its accessible porosity (to a methane molecule) can increase from 0 to ∼30%. We also show that these deformations are significantly nonaffine (i.e., nonhomogeneous), especially upon fluid adsorption or desorption.

Impact of Nanoporosity on Hydrocarbon Transport in Shales' Organic Matter.

In a context of growing attention for shale gas, the precise impact of organic matter (kerogen) on hydrocarbon recovery from unconventional reservoirs still has to be assessed. Kerogen's microstructure is characterized by a very disordered pore network that greatly affects hydrocarbon transport. The specific structure and texture of this organic matter at the nanoscale is highly dependent on its origin. In this study, by the use of statistical physics and molecular dynamics, we shed some new lights on hydrocarbon transport through realistic molecular models of kerogen at different level of maturity [ Bousige et al. Nat. Mater. 2016 , 15 , 576 ]. Despite the apparent complexity, severe confinement effects controlled by the porosity of the various kerogens allow linear alkanes (from methane to dodecane) transport to be studied only via the self-diffusion coefficients of the species. The decrease of the transport coefficients with the amount of adsorbed fluid can be described by a free volume theory. Ultimately, the transport coefficients of hydrocarbons can be expressed simply as a function of the porosity (volume fraction of void) of the microstructure, thus paving the way for shale gas recovery predictions.

From cellulose to kerogen: molecular simulation of a geological process.

The process by which organic matter decomposes deep underground to form petroleum and its underlying kerogen matrix has so far remained a no man's land to theoreticians, largely because of the geological (Myears) timescale associated with the process. Using reactive molecular dynamics and an accelerated simulation framework, the replica exchange molecular dynamics method, we simulate the full transformation of cellulose into kerogen and its associated fluid phase under prevailing geological conditions. We observe in sequence the fragmentation of the cellulose crystal and production of water, the development of an unsaturated aliphatic macromolecular phase and its aromatization. The composition of the solid residue along the maturation pathway strictly follows what is observed for natural type III kerogen and for artificially matured samples under confined conditions. After expulsion of the fluid phase, the obtained microporous kerogen possesses the structure, texture, density, porosity and stiffness observed for mature type III kerogen and a microporous carbon obtained by saccharose pyrolysis at low temperature. As expected for this variety of precursor, the main resulting hydrocarbon is methane. The present work thus demonstrates that molecular simulations can now be used to assess, almost quantitatively, such complex chemical processes as petrogenesis in fossil reservoirs and, more generally, the possible conversion of any natural product into bio-sourced materials and/or fuel.

Free Volume Theory of Hydrocarbon Mixture Transport in Nanoporous Materials.

Despite recent focus on shale gas, hydrocarbon recovery from the ultraconfining and disordered porosity of organic matter in shales (kerogen) remains poorly understood. Key aspects such as the breakdown of hydrodynamics at the nanoscale and strong adsorption effects lead to unexplained non-Darcy behaviors. Here, molecular dynamics and statistical mechanics are used to elucidate hydrocarbon mixture transport through a realistic molecular model of kerogen [ Bousige, C.; et al. Nat. Mater. 2016 , 15 , 576 ]. Owing to strong adsorption effects, velocity cross-correlations between the mixture components and between molecules of the same species are shown to be negligible. This allows estimation of each component permeance from its self-diffusivity, which can be obtained from single-component data. These permeances are found to scale with the reciprocal of the alkane length and decrease with the number of adsorbed molecules following a simple free volume theory, therefore allowing mixture transport prediction as a function of the amount of trapped fluid.

Pore network model of electrokinetic transport through charged porous media.

We introduce a method for the numerical determination of the steady-state response of complex charged porous media to pressure, salt concentration, and electric potential gradients. The macroscopic fluxes of solvent, salt, and charge are computed within the framework of the Pore Network Model (PNM), which describes the pore structure of the samples as networks of pores connected to each other by channels. The PNM approach is used to capture the couplings between solvent and ionic flows which arise from the charge of the solid surfaces. For the microscopic transport coefficients on the channel scale, we take a simple analytical form obtained previously by solving the Poisson-Nernst-Planck and Stokes equations in a cylindrical channel. These transport coefficients are upscaled for a given network by imposing conservation laws for each pores, in the presence of macroscopic gradients across the sample. The complex pore structure of the material is captured by the distribution of channel diameters. We investigate the combined effects of this complex geometry, the surface charge, and the salt concentration on the macroscopic transport coefficients. The upscaled numerical model preserves the Onsager relations between the latter, as expected. The calculated macroscopic coefficients behave qualitatively as their microscopic counterparts, except for the permeability and the electro-osmotic coupling coefficient when the electrokinetic effects are strong. Quantitatively, the electrokinetic couplings increase the difference between the macroscopic coefficients and the corresponding ones for a single channel of average diameter.

Numerical homogenization of electrokinetic equations in porous media using lattice-Boltzmann simulations.

We report the calculation of all the transfer coefficients which couple the solvent and ionic fluxes through a charged pore under the effect of pressure, electrostatic potential, and concentration gradients. We use a combination of analytical calculations at the Poisson-Nernst-Planck and Navier-Stokes levels of description and mesoscopic lattice simulations based on kinetic theory. In the absence of added salt, i.e., when the only ions present in the fluid are the counterions compensating the charge of the surface, exact analytical expressions for the fluxes in cylindrical pores allow us to validate a new lattice-Boltzmann electrokinetics (LBE) scheme which accounts for the osmotic contribution to the transport of all species. The influence of simulation parameters on the numerical accuracy is thoroughly investigated. In the presence of an added salt, we assess the range of validity of approximate expressions of the fluxes computed from the linearized Poisson-Boltzmann equation by a systematic comparison with LBE simulations.